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In this paper we construct integrable three-dimensional quantum-mechanical systems with mag- 
netic fields, admitting pairs of commuting second-order integrals of motion. The case of Carte- 
sian coordinates is considered. Most of the systems obtained are new and not related to the 
separation of variables in the corresponding Schrodinger equation. 



Let us consider the stationary Schrodinger equation for particle moving in external electromag- 
netic field in three-dimensional Euclidean space 



m 

Hip = Eip, H = -p 2 + V(x,y,z) + A i (x,y,z)pi+p i Ai(x,y,z), (1) 

■ 2 

c3 . 

where V(x,y,z) and A = (Ai{x, y, z), A^ix, y, z), A^[x, y, z)) are scalar and vector potentials 
of electromagnetic field, respectively. Here and below we use the notation p = —ih\7 and the 
summation from 1 to 3 over the repeated indices is understood. 

By analogy with classical Hamiltonian mechanics this system is called integrable if there 
CN \ exists a pair of quantum-mechanical operators P and Q which commute with each other as well 

as with the Hamiltonian H, i.e., the following relations hold 

<N 
(N 



[H,Q] = [H,P] = [P,Q]=0. 



Moreover, all three operators H, Q and P are algebraically independent, i.e., any of them cannot 
be represented as polynomial of two others [1]. 

In this paper we restrict ourselves by the case of Q and P being quadratic polynomials of p 

Q = ttik(x,y,z)pip k + f i (x,y,z)p i + j 1 (x,y,z), 
P = f3 ik {x,y,z)pip k + gi(x,y,z)pi + -f 2 {x,y,z), 

In classical mechanics integrable systems are interesting since their motion in a phase space 
is more ordered, namely, it is restricted to a torus. In quantum mechanics integr ability of an 
n-dimensional quantum-mechanical system, i.e., existing of n quantum integrals of motion (op- 
erators, which commute with each other as well as with the operator of equation ([I])), simplifies 
the problem of determining of energy spectrum and wave functions even when integrability does 
not lead to the variables separation, but only to the so-called "quasi-separation of variables" [2]. 

Therefore, the classification problem of all potentials V and A of electromagnetic field, for 
which quantum-mechanical problem is integrable in the sense described above, is of current 
importance. 

In the three-dimensional case of scalar potential, when there is no any magnetic field, this 
problem was solved as early as in 1967 by Yakov A. Smorodinsky with coauthors [3,4]. They 
have proved that if there exist two quantum integrals of motion of first or second order with 
respect to p (in the sense described above) then there exists a possibility of variable separation 
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in the corresponding Schrodinger equation [4]. The inverse statement is also true [3]. Thus, 
eleven classes of integrable potentials V obtained by them coincide with the results of classic 
paper of Eisenhart [5], where he as early as in 1948 has described all scalar potentials for which 
corresponding Schrodinger equation (or Hamilton-Jacobi equation in classical mechanics) admits 
variable separation at least in one of eleven coordinate systems. 

The next step was done in 1972 by V.N. Shapovalov with coauthors [6]. They have ob- 
tained complete classification list of vector-potentials with nonzero magnetic field, for which 
corresponding Schrodinger equation ([T]) admits variable separation at least in one of eleven co- 
ordinate systems, and list of corresponding pairs of operators which commute with each other 
and with operator of the equation. 

Posterior results in this direction are connected with works of P. Winternitz and his coau- 
thors [7-9] . For two-dimensional case they have discovered that existing of second-order integrals 
of motion in magnetic field does not guarantee possibility of variable separation. Nevertheless, 
even in this case, integrals of motion are classified on equivalence classes under action of Euclid- 
ian group and second-order terms with respect to pi in these integrals have the same form as 
in the case of pure scalar potential. It was also shown that in magnetic field quantum case [9] 
does not obviously coincide with classical one [7,8], namely, constructed vector-potentials can 
depend on Plank constant h in a nontrivial way. 

In this paper we make the next step in classification of potentials of electromagnetic field V 
and A in three-dimensional Euclidian space, for which corresponding quantum-mechanical sys- 
tem described by equation ([T]) is integrable in the sense explained above, i.e., for which there 
exists a pair of operators ([2]) that commute with each other as well as with the operator of 
the equation. As a result, we obtain a number of vector-potentials for which corresponding 
Schrodinger equation is integrable but does not admit variable separation. Therefore, these 
potentials did not appear in the Shapovalov's classification [6] and are new. 

At first we consider single operator Q of the form ([2]) which commutes with the operator H 
of Schrodinger equation dl|). The commutator [Q,H] contains terms of zero, first, second and 
third orders with respect to pi, which coefficients have to vanish. Coefficients of the third power 
of pi give the following system: 

da ik 

-K PmPiPk = 0. 

ox m 

After solving this system we obtain that the operator Q can be presented as symmetric bilinear 
polynomial of the infinitesimal generators of a group of motions of three-dimensional Euclidean 
space E%, i.e., symmetry group of the three-dimensional Schrodinger equation for free particle 
(the Helmholtz equation): 

Q = a ik MiM k + b ik (piM k + M k pi) + CijPip k + fi(x, y, z)pi + 71 (x, y, z), (3) 

where an-, bi k and Cj,- are constants, Mi is operator of rotation, namely, Mi = £ikl x kPh where 
£i k i is completely antisymmetric tensor. 

Comparing constructed form (j3J) of operator with similar form of operator in the case of pure 
scalar potential [4], we conclude that second-order terms with respect to pi remain the same 
after appearing of nonzero magnetic field (as well as in two-dimensional case [9]). 

Therefore, analogously to [4], a pair of commuting operators P and Q having the form ([3|) 
can be reduced by rotations and translations of coordinate system and by the transformations 

Q' = fiP + V Q + XH 

to one of eleven classes corresponding to eleven classical coordinate systems which provide 
separation of variable in the three-dimensional Schrodinger equation for free particle. 
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That is why in our case of nonzero magnetic field we obtain that the second-order terms with 
respect to pi in these eleven commuting pairs of the operators P and Q are of the same form as 
ones presented in [4]. However, in contrast to the pure scalar case, some or all of the coefficients 
of the first power of pi are nonzero functions. The forms of these functions determine the form 
of magnetic field (in scalar case they equal zero). 

In this paper we completely solve the simplest "Cartesian" case, i.e. the case 

Q = Pi + f{x,y,z)p + -yi(x, y,z), (4) 

P = pI +9(x,y,z)p + -f 2 (x,y,z). (5) 

Splitting the equations [H, Q] = [H, P] = [P,Q] = with respect to different powers of pi, 
we obtain an overdetermining system of PDEs on unknown functions fx, f 2 , f 3 , gi, g 2 , 93, li, 
72 and V, A u A 2 , A 3 . 

The coefficients of the highest powers of pi give the following system 

h = h{x), f 3 = f 3 (x), gx = gi{y), 93 = 93(9), hy = 92x, 
& 0) + fly = AA 2x , 93 (y ) + 92z = 4A 3y , AA lx = f lx , 
f 3 (x) + hz = 4^3x, 9i(y) + 92x = 4A ly , 4A 2y = g 2y . 

Its general solution for A is 

4,A X = s x + h x + g 1 (y) + n(z), 
AA 2 = s y + k 2y + f 2 (x) + r 2 (z), 
AA 3 = s z + k u + hz + / 3 (a;) + 93{y) + r' 3 (z), 

where s = s(x, y, z), k\ = k\(x, z) and k 2 = k 2 (y, z). 
The gauge transformation 

A^A + VF, F = s(x,y,z) + k 1 {x,z) + k 2 (y,z)+r 3 {z), 
simplifies the obtained expression for A in the following way 

Ai = \{gx(y) + ri(z)), A 2 = \{f 2 (x)+r 2 (z)), A 3 = ~(/ 3 (x) + g 3 (y)). 

With this expression for A in hand, we obtain from the coefficients of the lowest powers of 
Pi the system of ODEs on the functions g%(y), r\(z) f 2 (x), r 2 (z), f 3 (x) and g 3 (y): 

f2(x)g' 3 (y) = gi(y)f 3 (x), 

ri (z)f^x) = f 3 (x)r' 2 (z), (6) 
r 2 (z)g[(y) = gz{y)r' x (z). 

It is obvious that the Schrodinger equation with vector-potential (TjQ) is invariant with respect 
to permutations of A±, A 2 and A 3 , which are done simultaneously with permutations of the 
variables x±, x 2 and x 3 . Equations ([6]) are invariant with respect to permutations of the functions 
9i{y)i r i{ z )i h(x), r 2 (z), f 3 (x) and g 3 (y). These equivalence transformations can be represented 
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in such a way 



/ h 93 n \ 

h 9i r 2 

A 1 A 2 A 3 

\ x y z ) 

( h r 2 g\ \ 

/2 n 93 

At A 3 A 2 

\ x z y J 



( n f 2 g 3 \ 

r2 f3 91 

A 3 At A 2 

\ z x y ) 

( 9i h r 2 \ 

93 /2 rt 

A 2 At A 3 

\ y x z ) 



( 93 rt f 2 \ 

91 r 2 f 3 

A 2 A 3 At 

V y z x J 

( r 2 gt f 3 \ 

rt 53 h 

A 3 A 2 At 

\ z y x ) 



Usage of these equivalence transformations allows us to describe exhaustively all solutions of 
equations ([6]). This gives the complete description of all possible forms of vector-potentials A. 
The residual coefficients of powers of pt serve for determining of scalar component V of vector- 
potential and put additional constraints on the functions gt(y), rt(z), f 2 (x), r 2 (z), f 3 (x) and 

93 (y)- 

Below we adduce the final results of our calculations. Here and below denotes magnetic 
field, namely, 0, = rot A. 
Case 1. 

1=0, Q = 0, V = ut(x)+u 2 (y) + u 3 (z), Q = pl + 2ut(x), P = p 2 2 + 2u 2 (y). 

This case corresponds to zero magnetic field and is contained in the Eisenhart's classification [5]. 
According to his results all scalar potentials for which corresponding Schrodinger equations admit 
variable separation in Cartesian coordinates are exhausted by ones having the form 

V = ut(x) + u 2 {y) + u 3 {z). 

Case 2. 



/ vt(z) } 




f -w) \ 


\ V2(Z) 


v[(z) 1 









V = v 3 (z), Q 



P 



Case 3. 





, n= | 


' g'(y) \ 




: ) 


-f'(x) J 


, V = ut(x) + u 2 {y) 


f(x) + g(y) ) 









Q = pi + 4/(x)p 3 + 2ut(x), p = p 2 2 + Mv)p3 + 2 Mv)- 



The cases 2 and 3 where obtained by Shapovalov et al. [6]. According to their results these 
two cases exhaust all vector-potentials with nonzero magnetic field for which corresponding 
Schrodinger equations admit variable separation in Cartesian coordinates. 

Next three cases are not connected with variable separation and, therefore, these potentials 
did not appear in the Shapovalov's classification [6] and they are new. 



4 



Case 4. 

/ g'(y) \ ( o 

1= f{x) , n=[ o 

V / \f"(x)-g"(y) 

V = -(C 3 f(x) + C 3 g(y) + 2C 2 f 2 (x) + 2C l9 2 (y) + r(z) + Ag{y)f"{x) + 4/(x)y"(y)), 
Q=P 2 i+ 4/'(x)p2 - 2(4g(y)f"{x) + 2C 2 /(^) 2 + C 3 /(*)), 

P = p 2 + 4y'(y) Pl - 2{Af{x)g'\y) + 2dy(y) 2 + C 3 y(y)), 

where the functions f{x) and g(y) are solutions of the ODEs 

f(x) = Cf(x) + C l f(x) + C i , (7) 
y"(y) = Cy 2 (y) + C 2 y(y) + C 5 . (8) 

If C = these are linear second-order ODEs. The case C ^ is more interesting. If additionally 
Ci 7^ (resp. C 2 7^ 0), then solution of equation ([7]) (resp. ([8])) is first Painleve transcendent. If 
C = C\ = (resp. C = C 2 = 0), then ([7]) (resp. flS}) is the Weierstrass equation which solutions 
are expressed either via Weierstrass functions or via elementary ones depending on values of the 
parameter C4 (resp. C5) and integration constants. See, e.g. [10] for details. 
Case 5. 

/ g'{y) \ ( Cg'(y) 

A=[ f(x) , n=\ -Cf'(x) 

\ Cf(x) + Cg(y) J \ f"(x)-g"(y) 

V = -(C 3 f(x) + C 3 g(y) + 2C 2 / 2 (x) + 2dy 2 (y) + Ag(y)f"(x) + Af(x)g" (y)), 
Q=pl + 4(/ / (x)p 2 + Cf(x) P3 ) - 2(4g(y)f"(x) + 2C 2 /(x) 2 + C 3 f(x)), 
P=p 2 2 + 4(g'(y) P i + Cg{y)p 3 ) - 2(4f(x)g"(y) + 2C ig (y) 2 + C 3 g(y)), 

where the functions f{x) and g{y) are solutions of the ordinary differential equations 

f"( x ) = C 6 f(x) + C 1 f(x) + C 4 , 
g"(y) = C 6 g 2 (y) + C 2 g(y) + C 5 , 

which are integrated analogously to the previous case. 
Case 6. 



/ w'M + v'M \ ( w'i(y)-v' 2 \z) 

' u' 2 (x) + v' 2 (z) , n = -\ v'{{z) - u>'{x) 
u' 3 (x) + w' 3 (y) J \ u 2 (x) - w'{(y) 



V = --(iii(x) +w 2 (y) +v 3 (z) + w 1 (y)u 2 (x) + v 1 (z)u 3 (x) + 
u 2 (x)w'{{y) + v 2 (z)w 3 (y) + u 3 (x)v'{(z) + w 3 {y)v 2 (z)), 

Q = p\ + u' 2 {x)p 2 + u' 3 (x)p 3 - ^(wi(y)u 2 (x) + vi(z)u 3 (x) + Ui(x)), 
P = p\ + + w' 3 {y)p 3 - -ziu 2 {x)w'l{y) + ^(^^(y) + w 2 (y)), 
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where the functions u 2 (x), u 3 (x), wi(y), w 3 (y), v±(z), v 2 (z), u±(x), W2(y) and v%{z) are defined 
in a special way and described by the following four cases: 
Case 6.1. 

u 2 (x) = a 3 (ri cosh(aix) + k\ sinh(aix)), u 3 (x) = a 2 {r\ sinh(aix) + k\ cosh(aix)), 
w\{y) = a 3 (r 2 cosh(a 2 y) + k 2 sinh(a 2 y)), w 3 (y) = a\(r 2 sinh(a 2 y) + k 2 cosh(a 2 y)), 
vi(z) = a 2 (r 3 cosh(a 3 z) + k 3 sinh(a 3 z)), v 2 (z) = ai(r 3 sinh(a 3 z) + k 3 cosh(a 3 z)), 
a 2 a 2 

u\(x) = 2 ^ 3 ((r 2 + kf) cosh(2aix) + 2rik± sinh(2aix)) + C(r\ cosh(aix) + k\ sinh(aix)), 
a 2 a 2 

W2{y) = -^p ((r 2 + k 2 ) cosh(2a 2 y) + 2r 2 k 2 sinh(2a 2 y)) + C(r 2 cosh(a 2 y) + k 2 sinh(a 2 y)), 
a 2 a 2 

vs(z) = ((rf + kf) cosh(2a 3 ^) + 2r 3 k 3 sinh(2a 3 z)) + Ci(r 3 cosh(a 3 z) + k 3 sinh(a 3 z)), 

with 5 possible subcases: 

a) C = 0, Ci = 0; 

b) n = ki, r 2 = k 2 , r 3 = fc 3 , C\ = C; 

c) n = ki, r 2 = -k 2 , r 3 = -k 3 , C\ = C; 

d) n = -ki, r 2 = k 2 , r 3 = -k 3 , C\ = -C; 

e) n = -ki, r 2 = -k 2 , r 3 = k 3l C\ = -C. 

Case 6.2. 

u 2 (x) = a 3 (r± sin(aix) — k\ cos(aix)), u 3 (x) = a 2 (r\ cos{a\x) + k\ sin(aix)), 
w\{y) = a 3 (r 2 sin(a 2 y) - k 2 cos(a 2 y)), w 3 (y) = ai(r 2 cos(a 2 y) + k 2 sin(a 2 y)), 

= a 2 (r 3 cosh(a 3 z) + fc 3 sinh(a 3 z)), v 2 (z) = a\(r 3 sinh(a 3 z) + k 3 cosh(a 3 z)), 
a 2 a 2 

ni(x) = 2 ^ 3 ((r 2 — A; 2 ) cos(2aix) + 2ri&i sin(2aix)) + C{r\ sin(aix) — k\ cos(aix)), 
a 2 a 2 

i«2 (y) = -^p {(r 2 - k 2 ) cos(2a 2 y) + 2r 2 k 2 sin(2a 2 y)) + C(r 2 sin(a 2 y) - k 2 cos(a 2 y)), 
2 2 

1*3(2) = _^^2 ^ r 2 + fc 2) cosh (2a 3 z) + 2r 3 /c 3 sinh(2a 3 z)) + 
Ci(r 3 cosh(a 3 z) + /c 3 sinh(a 3 z)), 

with 5 possible subcases: 



a) 


C 


= 0, 


Ci 


= 0; 










b) 


n 


= iki, 


T2 


= ~ik 2 , 


r 3 


= ~k 3 , 


Ci 


= iC- 


c) 


ri 


= ih, 


T2 


= ik 2 , 


r 3 


= h, 


Ci 


= iC- 


d) 


n 


= -iki, 


T2 


= ~ik 2 , 


r 3 


= h, 


Ci 


= -iC; 


e) 


n 


= -iki, 


T2 


= ik 2 , 


r 3 


= -h, 


Ci 


= -iC. 
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Case 6.3. 



u 2 (x) = 03(^1 cos(aix) + k\ sin(aix)), u 3 (x) = ia 2 (ri sin(aix) — k\ cos(aix)), 
wi{y) = a 3 {r 2 cos(a 2 y) + k 2 sin(a 2 y)), w 3 (y) = iai(r 2 sm(a 2 y) - k 2 cos(a 2 y)), 
vi(z) = a 2 (r 3 cos(a 3 z) + k 3 sin(a 3 z)), v 2 (z) = iai(r 3 sin(a 3 z) - fo cos(a 3 z)), 



ui(x) = - 


n 2 n 2 

- f ((r? 


-kf) 


cos(2aix) + 2r\k\ sin(2aix)) + C(r\ cos(aix) + k\ sin(aix)), 


w 2 (y) = - 


_a 2 a§ ((r| 


~ k 2 ) 


1 cos(2a 2 y) + 2r 2 k 2 sin(2a 2 y)) + C(r 2 cos(a 2 y) + k 2 sin(a 2 y)), 


v 3 {z) = - 


n 2 « 2 

«1«2 // 2 

4 ^ r 3 " 


-kl) 


cos(2a 3 z) + 2r 3 k 3 sin(2a 3 ,z)) + C\(r 3 sin(a 3 z) — k 3 cos(a 3 z)), 


with 5 possible 


subcases: 







a) 


C 


= 0, 


Ci = 


:0; 










b) 


ri 


= ifo, 


r 2 = 


-ifo, 


r 3 


= ifo, 


Ci 


= iC; 


c) 


ri 


= -ih, 


r 2 = 


-ik 2 , 


r 3 


= -ifo, 


Ci 


= iC; 


d) 


n 


= -ifo, 


r 2 = 


ik 2 , 


r 3 


= ik 3 , 


Ci 


= -iC; 


c) 


ri 


= ifo, 


r 2 = 


ik 2 , 


r 3 


= -ifo, 


Ci 


= -iC. 



Case 6.4. 

u 2 (x) = a 3 (ri cosh(aix) + k\ sinh(aix)), u 3 (x) = — ia 2 {r\ sinh(aix) + k\ cosh(aix)), 
wi{y) = a 3 {r 2 cosh(a 2 y) + k 2 sinh(a 2 y)), w 3 (y) = -iai(r 2 sinh(a 2 y) + k 2 cosh(a 2 y)), 
v\(z) = a 2 (r 3 cos(a 3 z) + k 3 sin(a 3 z)), v 2 (z) = ia\(r 3 sin(a 3 ,z) - k 3 cos(a 3 z)), 

a 2 a 2 

u\(x) = 2 ^ 3 ((r 2 + kf) cosh(2aix) + 2r\k\ sinh(2aix)) + C(r\ cosh(aix) + k\ sinh(aix)), 
aV 

w 2 (y) = — — - ((r 2 + k\) cosh(2a 2 y) + 2r 2 A; 2 sinh(2a 2 y)) + C(r 2 cosh(a 2 y) + k 2 sinh(a 2 y)), 
a 2 a 2 

v 3 (z) = -1— ((r| - kf) cos(2a 3 z) + 2r 3 A; 3 sin(2a 3 z)) + Ci(r 3 sin(a 3 z) - k 3 cos(a 3 z)), 
with 5 possible subcases: 



a) 


C 


= 0, 


Ci = 


■0; 










b) 


r\ 


= ~h, 


r 2 = 


-k 2 , 


r 3 


= -ifo, 


Ci 


= C; 


c) 


ri 


= h, 


r 2 = 


-k 2 , 


r 3 


= ifo, 


Ci 


= C; 


d) 


ri 


= k\, 


r 2 = 


k 2 , 


r 3 


= -ifo, 


Ci 


= -C; 


e) 


n 


= ~ki, 


r 2 = 


k 2 , 


r 3 


= ifo, 


c\ 


= -C. 



Therefore, we obtain a number of new vector-potentials, for which the corresponding Schrodinger 
equation ([I]) is integrable in the sense described above. 

The author is grateful to Prof. P. Winternitz for useful discussions. This work was partly 
supported by a grant from NATO. 
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